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Strongly correlated electron systems with degenerate 
orbitals have attracted much interest. One of typical 
examples is the transition metal compound L1V2O4, 1 ' 
where the heavy-fermion behavior was observed down to 
very low temperatures. It was suggested that tig or- 
bitals as well as geometrical frustration are important 
to understand the ground state stabilized against other 
ordered states. 2 ' Another example is the /-electron sys- 
tem CeTIn 5 (T = Co, Ir, RI1), 3 - 5 ) where the supercon- 
ducting state competes with the magnetically ordered 
state. It was claimed that the effective Hamiltonian 
for this compound may be described by the degenerate 
Hubbard model. 6 ' Therefore, it is desirable to investi- 
gate low-energy properties of the degenerate Hubbard 
model systematically. In our previous paper, 7 ' we stud- 
ied the ground-state properties of the degenerate Hub- 
bard model by means of the dynamical mean field the- 
ory (DMFT) combined with the exact diagonalization 
method. 8 ' It was found that if the intra- and inter- 
band Coulomb interactions are nearly equal, the metallic 
ground state is stabilized up to fairly large interactions 
due to orbital fluctuations. Shown in Fig. 1 is a more 
precise phase diagram newly obtained by DMFT with 
the exact diagonalization. However, we were not able 




Fig. 1. Quasi-particle weight for the degenerate Hubbard model 
calculated by the DMFT within the exact diagonalization. Z = 
(Z ^ 0) corresponds to the insulating (metallic) phase. 

to systematically deal with dynamical properties in the 
previous paper because of the limitation inherent in the 
exact diagonalization of small clusters. We address this 
problem here by combining the quantum Monte Carlo 
simulations with DMFT. This also enables us to explore 
finite-temperature properties, which are useful to reveal 
the many-body character in correlated electron systems. 

We consider a correlated electron system with two- 
fold degenerate orbitals, which may be described by the 



Hubbard Hamiltonian, 



U ^niafUiai + U' ^2 n il<? n i2o> (1) 



where c iaiy (ci aa ) creates (annihilates) an electron with 
spin <x(=j,JJ and orbital a(— 1,2) at the ith site, and 
n iaa — c l aa c iaa- e fc represents the dispersion relation, 
fi the chemical potential, U and U' the intra-band and 
inter-band Coulomb interaction. We use the semicircu- 
lar density of states p(x) — 2/irD^Jl — (x/D) 2 for sim- 
plicity. We focus here on the symmetric case by setting 
fj, = hU + U' . In the following, we take the band width 
D as unity. 

To discuss how two kinds of Coulomb interactions af- 
fect the stability of the metallic state, we use DMFT, 8 ' 
which can incorporate local electron correlations pre- 
cisely. This method has been applied to the degener- 
ate Hubbard model by several groups. 7, 9 ~ 12 ' In this pa- 
per, we exploit the quantum Monte Carlo simulations 
based on the Hirsch-Fye algorithm 13 ' to clarify finite- 
temperature properties of the degenerate model. 

In Fig. 2, we show the density of states deduced by 
applying the maximum entropy method to the Monte 
Carlo data. When U' = 0.0, the system is reduced to the 
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Fig. 2. Density of states for the degenerate Hubbard model (U = 
3.0). The data are for inverse temperature T~ l = 1,2,4,8 and 
16 from the top to the bottom. 

single-band Hubbard model, where the metal-insulator 
transition occurs at a critical point U c — 2.94. 14 ' In 
the case {U,U') = (3.0,0.0), the system is in the in- 
sulating phase close to the transition point. Neverthe- 
less we can clearly observe the formation of the Hubbard 
gap with the decrease of temperature T. Introducing 
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the interband Coulomb interaction U', the sharp quasi- 
particle peak appears around the Fermi surface in the 
case (U,U r ) — (3.0,3.0). This implies that orbital fluc- 
tuations induced by the interband Coulomb interaction 
drives the system to the metallic phase. In fact, the lo- 
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Fig. 3. Local charge susceptibility as a function of the interband 
Coulomb interaction with a fixed U = 3.0. 

cal charge susceptibility is increased up to U ~ U' at low 
temperatures, as seen in Fig. 3. These results are con- 
sistent with the previous results. 7,9 ' 12 -* Further increase 
in the interband interaction suppresses spin fluctuations, 
leading the system to another type of the Mott insulator 
in the region of U' > U (see Fig. 1). 

To characterize the nature of two types of the Mott 
insulators, we investigate the temperature dependence 
of the local spin and orbital susceptibilities, which are 
defined as, 

Jo 

Xo = / dT<{ni(0)-n2(0)}{ni(r)-n2(T)}>^2) 
Jo 

where (3 = T _1 , n a = ]C Q ?W, n a = ^ CT n ai<7 and T 
is an imaginary time. We show the results obtained by 
the quantum Monte Carlo simulations within DMFT in 
Fig. 4. Let us first look at Figs. 4 (a) and (b) for 
U' = (equivalent to the single band model). Since the 
introduction of the intra-orbital interaction U sharpens 
the quasi-particle peak, this gives rise to the increase of 
the spin susceptibility Xs at low temperatures. On the 
other hand, the formation of the Hubbard gap suppresses 
not only the charge susceptibility, but also the orbital 
susceptibility \o- As seen from Figs. 4 (c) and (d), we 
encounter quite different behavior in the susceptibilities, 
when the interband interaction V is increased. Namely, 
the spin susceptibility is suppressed, while the orbital 
susceptibility is enhanced at low temperatures. We have 
checked that this tendency holds for larger U' beyond 
the condition U' ~ U. 

Therefore in the metallic phase close to the Mott insu- 
lator in the region of U > U'(U < U'), spin (orbital) fluc- 
tuations are enhanced whereas orbital (spin) fluctuations 
are suppressed with the decrease of temperature. These 
analyses shed some light on the reason why the metallic 
phase is particularly stable along the line U — U' . On 
this line, spin and orbital fluctuations arc almost equally 
enhanced, and this subtle balance is efficient to stabilize 
the metallic phase. When the ratio of interactions devi- 
ates from this condition, the system prefers either of two 




Fig. 4. Local spin and orbital susceptibilities as a function of tem- 
perature. 



types of the Mott insulating phases. 

In this paper, we have restricted our discussions to 
the paramagnetic phase. As mentioned above, spin 
and/or orbital fluctuations are particularly enhanced in 
the metallic phase close to the Mott insulator, which 
would trigger a quantum phase transition to a certain 
ordered phase. 10 ) This problem is to be addressed in the 
future work. Nevertheless, we think that the general fea- 
ture about the stability of metallic phase near U ~ V 
may hold even in those cases. 
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